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Radiant-Heat Transfer Between Nongray Parallel Plates 1 

Stanley Goodman 

Radiant-heat transfer rates between infinite parallel plates were computed for Inconel 
and aluminum over a wide range of plate temperatures. Results indicate that assuming 
the plates to be graybodies introduces an error of 2 to 29 percent in the computed heat- 
transfer rate. 

In cross-sectional heat flow through box beams typical of aircraft structures, it was 
found that radiant-heat transfer predominates in Inconel beams and is significant in 
aluminum beams. 



1. Introduction 

The current need for predicting thermal stresses 
and deformations in aircraft structures gives rise 
to a need for determining the transient temperature 
distribution in these structures. As the tempera- 
tures in a structure increase, radiation becomes 
more important as a mode of heat transfer. 

In computing radiant-heat transfer, a conventional 
assumption used is that the materials involved are 
graybodies, i. e., have spectral emissivities that are 
independent of wavelength. However, for polished 
metallic surfaces, spectral emissivities generally vary 
with wavelength, as is indicated by theoretical rela- 
tionships for spectral emissivity of smooth-surfaced 
electrical conductors fl] 2 and by published experi- 
mental data [2 to 7]. 

The purposes of this investigation were (1) to 
determine the magnitude of error caused by assuming 
that certain metals and alloys are graybodies in 
computing heat-transfer rates for a range of tem- 
peratures, temperature differences, and emissivities, 
and (2) to compare the radiant-heat transfer rate 
with the conductive-heat transfer rate in a box 
beam for a range of temperatures and temperature 
differences for two structural metals of different 
emissivities and thermal conductivities. Gaseous 
heat transfer is neglected. 



2. Symbols 



W- 



= Absorptivity 

= Emissivity 

= Thermal conductivity (Btu/hr ft °F) 

= Wavelength (microns) 

= Heat transfer rate (Btu/hr ft 2 ) 

= Heat transfer per unit time (Btu/hr) 

= Ratio of radiant-heat transfer rate to conductive 

heat-transfer rate 
= Electrical resistivity (ohm cm) 
= Stefan-Bcltzmann constant = (0. 1713 X10" s Btu/ft 2 

hr °R 4 ) 
-Temperature (°C or °F) 
= Temoerature (°R) 
-Micron (0.0001 cm) 
-Radiant flux density of a blackbody (Btu/hr ft 2 ) 



1 This work was conducted at the National Bureau of Standards under the 
sponsorship and with the financial assistance of the Office of Naval Research. 

2 Figures in brackets indicate the literature references at the end of this paper. 



Subscripts 

6 = Black body 
c= Conductive 

h — Hemisp herical 

X= Spectral 

= Normal 

r = Radiant 

t= Total 

1= Hotter of two parallel surfaces involved in radiant 

heat exchange 
2 = Cooler of two parallel surfaces involved in radiant 

heat exchange. 

3. Computations 

3.1. Materials 

Computations were limited to appropriate mate- 
rials for which the infrared spectral-emissivity curves 
are available in the literature. Such curves are 
presented in figure 1 for 99.6-percent-pure aluminum 
[7] and for Inconel [6]. 

3.2. Structure 

The structure analyzed was a box beam with two 
opposite walls at different uniform temperatures. 
The other two walls were assumed to be perfect 
reflectors with identical temperature distributions. 
Two beams were considered, one relatively thin- 
walled and one relatively thick-walled (see sketch, 
table 1). 

4. Method of Analysis 

For internal radiant-heat transfer the box beam 
described is obviously equivalent to a pair of infinite 
parallel planes with temperatures equal to the box- 
beam wall temperatures. 

In order to cover the entire wave band over which 
significant radiant-energy exchange occurs, values of 
e\ beyond the experimental wavelength ranges [6 
and 7] in figure 1 were obtained by extrapolation. 
The wave bands used in the computations were : for 
aluminum from \=0 to \=25 /z, and for Inconel 
from \=0 to X=15 jjl. For all cases considered, at 
least 95 percent of the emitted radiant energy is 
within these wavelength ranges. The proportion of 
this emitted radiant energy within the wave bands 
in which the experimental values of e were determined 
ranges roughly from 65 to 85 percent. 
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Figuke 1. Spectral normal emissivities of polished aluminum 
and as-rolled Inconel. 



To approximate the values of ex outside the range 
of experimental determinations, ex was extrapolated 
linearly for Inconel. For aluminum, e\ was taken 
as constant in the shorter wavelength range. (Rel- 
atively very little energy is emitted by the alum- 
inum in this wavelength range.) For the longer 
wavelengths, ex was assumed to vary inversely with 
Vx [1]. The extrapolated values are shown in 
figure 1. 

4.1. Computation of Radiant-Heat Transfer, 
Using Graybody Assumption 

In this case, the heat-transfer rate was computed 
from the formula given in reference [8] for net heat 
exchange between infinite parallel gray planes, 



2,=<r(T}-rt) 



1 



!/«!«»+ 1/«2«»— 1 



(1) 



To compute e th , e tQ was first computed from the 
measured values of e\ for aluminum and Inconel 
(fig. 1), by use of the relationship 



Table 1. Ratio, R, of radiant-heat transfer rate, q r , to conductive-heat transfer rate, q e , 
for thin-walled and thick-walled beams 



(VERTICAL WALLS ARE 
PERFECT REFLECTORS) 



— ■- t—c-ro 






• Total hemispherical emissivity. 



*» Thermal conductivity. 
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Material and material properties at temperature U 


Temperatures 


r R 


R 


Material 


& €th 


H 


eth 


h 


tl-t 2 


r=0.01 


r=0.03 








k 






b 


6=1 ft 


b=Hft 


1 


2 


3 


4 


5 


6 


7 


8 


9 






Btu/hrft°F 


hrft°F/Btu 


°C 


°c 


ft-i 








/ 0. 102 


133 


0. 767 X10 3 


300 












.115 


144 


.799 


400 


100 


0. 0019 


0.19 


0.021 




.130 


155 


.839 


500 


100 


.0033 


.33 


.037 


Aluminum 


( .130 
.113 


155 
166 


.839 
.681 


500 
600 


200 
100 


.0027 
.0054 


.27 

.54 


.030 
.060 






.113 


166 


.681 


600 


200 


.0040 


.40 


.044 




I .113 


166 


.681 


600 


300 


.0034 


.34 


.038 




( .445 
.504 


10.4 
10.8 


42.8 
46.7 


300 
400 


"""166 








.19 


19 


2.1 




.547 


11.5 


47.6 


500 


100 


.31 


31 


3.4 




.547 


11.5 


47.6 


500 


200 


.25 


25 


2.8 




.584 


12.1 


48.3 


600 


100 


.47 


47 


5.2 




.584 


12.1 


48.3 


600 


200 


.39 


39 


4.3 


Inconel _. _ _ 


{ .584 
.643 


12.1 
13.1 


48.3 
49.1 


600 
800 


300 
100 


.33 

.90 


33 

90 


3.7 
10.0 






.643 


13.1 


49.1 


800 


300 


.69 


69 


7.7 




.643 


13.1 


49.1 


800 


500 


.52 


52 


5.8 




.690 


14.1 


48.9 


1, 000 


100 


1.54 


154 


17.1 




.690 


14.1 


48.9 


1,000 


300 


1.24 


124 


13.8 




I .690 


14.1 


48.9 


1,000 


500 


.99 


99 


11.0 



e ' 0= l 



\ M ^-d(\T)-- 

-aT s 



X 



'' b\h 7/ A rji\ 
0~1 



(2) 



measure of the relative importance of radiant and 
conductive heat transfer in a box beam: 



Values of Wb\h/aT 5 SiS a function of XT 7 were taken 
from table 1 of reference [9] and substituted in eq 
(2), which was then integrated numerically. For 
aluminum, values of e th were obtained from e t0 
values by use of Schmidt and Eckcrt's relationship 
between e t0 and c^/e ro for smooth-surfaced metals [1]. 
Inconel (0.44<^ <C0.690), which forms an oxide 
surface layer on heating, was assumed to be inter- 
mediate between an electrical conductor and an 
insulator. The ratio e th je m should lie between 0.95 
and 1.05 [1]; e m was taken as equal to e m . Appro- 
priate values of e th were substituted in eq (1), which 
was then solved for g r . 

4.2. Computation of Radiant-Heat Transfer Rate 
From e\h 

Taking a Xh =e\ h and summing the scries of reflec- 
tions back and forth between planes 1 and 2 for the 
net radiant-heat transfer rate from plane 1 to plane 
2 for monochromatic radiation, g\, gives 



q\[(W bXh ) T = f-WWr = t 2 ] 



1 



1/61X^+1/62X^—1 



(3) 



For aluminum, e\ ?l was computed from Schmidt and 
Eckert's relationships between ex n and r e /\, and be- 
tween r e j\ and e\ h [1]. Values of W h \ h were com- 
puted from table 1 of [9]. For Inconel, because 
tth~eto, it follows that e\ h ~ e\ Q , and €\ h was taken as 
equal to ex . 

To compute the radiant-heat transfer rate be- 
tween planes 1 and 2, the appropriate values of 
W b \ h and e\ h were substituted in eq (3) and the equa- 
tion was then integrated numerically over the wave 
band X=0 to X=25 p, for aluminum and \=0 to 
X=15 m for Inconel. 

4.3. Computation of Conductive-Heat Transfer Rate 

Walls 1 and 2 of the box beam were assumed to 
have uniform temperature distributions. The heat 
transfer per hour was computed from 



Q - kA (kzIi), 



(4) 



where 

Q c = conductive-heat transfer per unit time 

(Btu/hr), 
A= total cross-sectional area of conducting 

walls (ft 2 ), 
b = distance between radiating wall surfaces 

(ft). 

Values of Jc, taken from reference [8] for aluminum 
and from reference [10] for Inconel, are listed in col- 
umn 3 of table 1 . The values of k used in eq (4) are 
taken at the temperature l/2(^ + ^ 2 )- 

Values of the ratio R^QrlQc^qlqc were used as a 



R-- 



Q r= ^ q r la q r b 

~Qc /j2 ^i-* 2 2^-^) r 



(5) 



where 



Z=length of beam (ft), 
a = width of radiating surface (ft), 
c=ra= thickness of each conducting wall (ft), 

b = distance between radiating wall surfaces 

(ft). 

5. Results 

Values of q r computed from the actual spectral 
emissivities, values of q r computed from total emis- 
sivities, using a graybody assumption, and the per- 
centage difference between corresponding values are 
listed in columns 4, 5, and 6, respectively, of table 2. 

The high temperatures for aluminum were in- 
tended to give an upper limit for the effect of radiant- 
heat transfer. 



Table 2. Computed radiant-heat flow rate, q r , between infinite 
parallel plates at temperatures t\ and t 2 



Material 


h 


t\-h 


a Qr 


h Qr 


a qr- b qr 












aq r 


1 


2 


3 


4 


5 


6 




°C 


°C 


Btu/hr ft* 


Btu/hr ft* 


% 




400 


100 


97 


95 


2 




500 


100 


178 


157 


12 


Aluminum 


500 
600 


200 

100 


283 
311 


271 
259 


4 

17 




600 


200 


451 


385 


15 




I 600 


300 


548 


483 


12 




/ 400 


100 


726 


543 


25 




500 


100 


1,250 


973 


22 




500 


200 


1,980 


1,460 


26 




600 


100 


1,980 


1,590 


20 




600 


200 


3, 240 


2, 510 


23 


InconeL-- 


600 


300 


3,970 


2,880 


27 




800 


100 


4, 140 


3, 540 


15 




800 


300 


9,110 


7,320 


24 




800 


500 


11,090 


7,820 


29 




1, 000 


100 


7, 660 


6, 780 


12 




1,000 


300 


17, 900 


15, 000 


16 




I 1, 000 


500 


22, 800 


17, 900 


21 



» Computed from spectral emissivity. & Computed from total emissivity. 

Ratios of radiant-heat transfer rate to conductive- 
heat transfer rate are listed in table 1 for a box beam 
in which all radiant-heat transfer is between the 
horizontal walls, and all conductive-heat transfer 
is in the vertical walls. Column 7 of the table lists 
ratios of radiant-heat transfer to conductive-heat 
transfer, R, times dimension ratio, r/b, where r=c/a 
is the ratio of vertical-wall thickness to horizontal- 
wall width, and b is the distance between horizontal- 
wall surfaces. Values of R for a box beam with 
c / a =0.01 and b=l ft are listed in column 8 of table 1, 
and values of R for a box beam with c/a=0.03 and 
6=1/3 ft are listed in column 9. 
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6. Discussion 

The shapes of the spectral-emissivity curves (see 
fig. 1) for aluminum and Inconel differ widely; for 
Inconel the curve decreases linearly with increasing 
wavelength. DeCorso and Coit [6] state that the 
spectral emissivity of Inconel should not vary with 
temperature if the surface condition is maintained. 
For aluminum the curve approaches graybody 
emissivity for \>2.5 \x, and is temperature dependent. 
The effect of t u (t 1 —t 2 ) } and e on the size of the error 
that is caused by assuming the materials are gray- 
bodies when computing radiant-heat flow rates is 
dependent on the shape of the spectral-emissivity 
curve. For aluminum the error increases with 
increasing t x and with decreasing {U — U). For 
Inconel, with greater e x (and therefore a\) than alumi- 
num, the error is greater than in aluminum and in- 
creases with decreasing ^ and with increasing (ti—t 2 ) . 

For aluminum and those materials that approach 
graybodies for the longer wavelengths, the graybody 
assumption may be used when the proportion of 
energy radiated at the shorter wavelengths is small. 
Pure aluminum appears to act like a graybody when 
its radiating surfaces are at temperatures lower than 
400° C. Due to the paucity of spectral-emissivity 
data, it is not known whether or not pure aluminum 
and Inconel have e x curves typical of other aircraft 
structural materials. 

The total emissivity of a material depends more 
on its surface condition, history of heating, oxidation, 
cleanliness of surface, and wetness than on its chemi- 
cal composition. For example, Hase [7] gives the 
following values of e x (and therefore a x ) for pure 
aluminum at 400° C for various surface conditions 
at \=2: Sandblasted, 0.63; dull, 0.29; polished, 0.17. 
For Inconel, DeCorso and Coit [6] give the following 
values of e t at 1,000° F after various periods of heat- 
ing: Initial, 0.28; after 15 min at 1,500° F, 0.53; 
after 30 min at 1,800° F, 0.62; after 135 min at 2,000° 
F, 0.80. It follows that for the cases considered, 
the error in predicting q r due to the error in the 
value of e may be greater than the error due to use 
of a graybody assumption for nongray materials. 

The beams chosen for analysis are common to 
many aircraft structures. Radiant-heat transfer 
will have its greatest value for this type of geometri- 
cal configuration and will be independent of the 
distance between the radiating walls. Because the 
quantity of energy radiated from a body varies as 
the fourth power of the absolute temperature of the 
body, whereas the quantity of energy conducted in a 
body depends on the temperature gradient, the ratio 
of radiant- to conductive-heat transfer rate, q r /q c , 
must be greatest for high temperatures. For box 
beams of the same geometry, q r /q c must increase 
with increasing e/k. The large differences between 
the q r /q c values for corresponding wall temperatures 
in aluminum and Inconel are due to the fact that e/k 
is comparatively small for aluminum and large for 
Inconel. With increasing temperature, e th /k gen- 
erally increases in the reasonable structure temper- 
ature range (see table 1). 

For a given structure, the size of the error in com- 



puted temperature distribution due to error in the 
computed radiant-heat transfer rate depends on the 
proportion of radiant-heat transfer to total heat 
transfer. It would be greater for Inconel, with a 
large e/k ratio, than for aluminum, with a small e/k ratio. 

7. Conclusions 

The error due to assuming pure aluminum and 
Inconel are graybodies in computing the radiant- 
heat transfer between infinite parallel plates of these 
materials over a wide range of temperatures and 
temperature differences varies from 2 to 29 percent. 
An error of this magnitude might also be encountered 
in determining the emissivity of a given material, 
because of the large variation in e with surface con- 
ditions. The size of the error for given plate tem- 
peratures and emissivities depends on the shape of 
the spectral-emissivity curve. Lack of spectral- 
emissivity data makes it difficult to determine the 
error for other aircraft materials. 

There is a need for determination of spectral 
hemispherical emissivities of aircraft materials at 
various temperatures and for various surface con- 
ditions. 

The ratio of radiant-heat transfer to conductive- 
heat transfer in a pure aluminum box beam (in the 
absence of gaseous heat transfer) ranges from 
0.19 to 0.54 for a thin- walled beam and from 0.021 
to 0.06 for a thick-walled beam over a range of hot- 
wall temperatures from 400° to 600° C and a range of 
wall-temperature differences from 100° to 300° C. 
For an Inconel box beam the corresponding ratios 
are: thin-walled, 19 and 154; thick-walled, 2.10 and 
17.1, over a range of hot-wall temperatures from 
400° to 1,000° C and a range of wall-temperature 
differences of 100° to 500° C. These results for a 
particular box beam indicate that radiant-heat 
transfer predominates for Inconel structures and is 
significant for thin-walled aluminum structures, at 
elevated temperatures. 
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